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Abstract: 

The mode-fluctuation distribution PiW) is studied for chaotic as well as for non-chaotic 
quantum billiards. This statistic is discussed in the broader framework of the E(k, L) functions 
being the probability of finding k energy levels in a randomly chosen interval of length L, and the 
distribution of n(L), where n(L) is the number of levels in such an interval, and their cumulants 
Ck(L). It is demonstrated that the cumulants provide a possible measure for the distinction 
between chaotic and non-chaotic systems. The vanishing of the normalized cumulants Cjt, 
k > 3, implies a Gaussian behaviour of PiW), which is realized in the case of chaotic systems, 
whereas non-chaotic systems display non-vanishing values for these cumulants leading to a 
non-Gaussian behaviour of PiW). For some integrable systems there exist rigorous proofs of 
the non-Gaussian behaviour which are also discussed. Our numerical results and the rigorous 
results for integrable systems suggest that a clear fingerprint of chaotic systems is provided by 
a Gaussian distribution of the mode-fluctuation distribution PiW). 
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I Introduction 



One of the main research lines of quantum chaology is devoted to the search of fingerprints left 
on a quantum system by its classically chaotic counterpart. In classically bounded conservative 
systems, which are chaotic, one of the most obvious properties of chaos manifests itself in the 
long-time behaviour, like the exponential decay of certain correlation functions for t — > oo. Their 
quantum mechanical counterparts do not display such a behaviour since the time-evolution 
operator U(t) = e~fr Ht , H being the Hamiltonian, leads to an almost periodic correlation 
function (cf)\U(t)(f)') = J2^=i dne~^ Ent , where </>, <// denote two states and E n the quantal levels of 
H. Thus, with respect to the long-time behaviour, there does not exist a fingerprint reflecting 
the classical properties. Therefore, in quantum chaology the limit t = oo is studied instead of 
the long-time asymptotics t — > oo, which leads to the study of eigenstates, i.e., wave functions 
and quantal energy spectra, especially of their statistical properties. 

This paper is devoted to a comparison between the statistics of the quantal levels of integrable 
and chaotic systems. For integrable systems, the short-range correlations in quantal energy 
spectra are, in general, numerically found to be Poissonian, but the long-range correlations 
show saturation effects [1, 2, 3]. In contrast, in chaotic systems the short-range correlations are 
numerically found to be in accordance with random-matrix theory (RMT) [4, 5, 6, 7], whereas 
long-range correlations experience modifications leading again to saturation effects. This general 
behaviour of chaotic systems is violated by the class of arithmetical systems [8, 9, 10, 11, 12, 13], 
where, e.g., the level-spacing distribution and the two-point statistics nearly behave as for 
classically integrable systems. 

In most numerical studies concerning the statistical properties of quantal levels of chaotic 
systems, the nearest-neighbour level-spacing distribution P(0,L), the number variance S 2 (L) 
or the spectral rigidity A 3 (L) is studied. However, in this paper we will focus on the E(k, L) 
functions, which give the probability to find k levels in a randomly chosen interval of length 
L. They constitute a fundamental measure, since other statistics like the three just mentioned 
above are completely determined by them. Furthermore, their numerical calculation is simple 
and well defined, since no binning is necessary. In addition, as a function of L one obtains 
smooth curves and significant results even in cases where only a small number of eigenvalues is 
known. The E(k, L) functions are related to the moments and cumulants of the distribution of 
n(L), which in turn are important for the behaviour of the mode-fluctuation distribution PiW) 
[14, 15]. It now seems that the latter provides a clear fingerprint of quantum chaos. 

The paper is organized as follows. In section II we provide some general remarks on the 
eigenvalue statistics of billiard systems, since we restrict our numerical analysis to those systems. 
Section III gives a summary of important properties of the E(k, L) functions as well as of their 
behaviour in random matrix theory, and some phenomenological properties. Section IV discusses 
n(L) and its corresponding distribution. The distribution of the normalized mode-fluctuations 
W(x) and its connection to the distribution of t]l(x) is discussed in section V. Exact results 
for mode-fluctuation distribution PiW) in the case of integrable systems are summarized in 
section VI. Finally, section VII gives the numerical comparison between non-chaotic and chaotic 
systems. 
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II General remarks on eigenvalue statistics in the case of 
quantum billiards 



The systems to be studied are either given by the free motion of a point particle inside a compact 
Euclidean domain with elastic reflections at the boundary <9fi, or by the geodesic motion on a 
compact Riemann surface with constant negative curvature, or by billiards on compact domains 
with constant negative curvature. In each case the quantum mechanical system is determined 
by the stationary Schrodinger equation (h = 2m = 1) 

- A* n (q) = EA(q) , (I) 

with the appropriate Dirichlet, Neumann or periodic boundary conditions, respectively. A 
denotes the Laplace-Beltrami operator, which on the Euclidean plane is the usual Laplacian. 
Billiard systems belong to the class of scaling systems, which means that one has the scaling 
property S^(XE) = \ a S~ / (E) of the action S^(E) = J p c/q of periodic orbits 7. If one defines 
x = E a } one has S^(E) = (E / E ) a S 1 (E ) =: xRy where i? 7 is independent of E. Also geodesic 

flows on Riemannian manifolds are scaling systems, where x = p = \J E — 1/4 is the modulus 
of the momentum. Hamiltonian systems with potentials scaling according to V(Aq) = A K V(q), 
A > 0, provide another class of scaling systems, where a = | + K For billiards the semiclassical 
limit corresponds to the high-energy limit E n — > 00. 

Due to the compactness of 0, the quantal energy spectrum {E n } is purely discrete; thus one 
can define the spectral staircase function N(E) (mode number, integrated level density) 

N(E) := #{n\E n < E} , (2) 

which counts the number of energy levels E n below a given energy E. N(E) can be separated 
into a mean smooth part N(E) and a fluctuating part, the mode fluctuations Nr uc (E) 

N(E) = 77(E) + N nuc (E) . (3) 

For two-dimensional billiards, N(E) is given by the generalized Weyl's law [16] and an additional 
term N\ ons (E) 

N(E) = ^E- ^VE + C + N loag (E) , (4) 

where A denotes the area of the billiard, and C := C~ — C + , where C~ and C + are the lengths 
of the boundary dQ with Dirichlet and Neumann boundary conditions, respectively. C takes 
curvature and corner corrections into account. The term N\ ong (E) is zero or negligible for most 
systems, but for some systems N(E) possesses a slowly varying long-range oscillation which we 
denote by N\ ong (E). This was first observed [17] in the case of the stadium billiard [18, 19] (see 
sec. VII. 2), where by subtracting the contribution of the so-called bouncing-ball orbits [20] the 
expected behaviour of the spectral statistics could be restored. Further contributions to N\ ong (E) 
in the case of the stadium billiard are higher order corrections [20, 21] caused by edge-orbits, 
which run into the point where circle and straight line join. The same applies to the truncated 
hyperbola billiard (see section VII. 2), where a prominent contribution to N(E) is given by 
families of closed non-periodic orbits running into a boundary point where the curvature is 
discontinuous [22]; but here their contribution vanishes for E — > 00, whereas the contribution of 
the bouncing-ball orbits in the stadium billiard increases in an unbounded manner. The main 
difference between N\ ong (E) and A^fl uc (£^) is that the oscillations in N\ ong (E) are regular whereas 
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the fluctuations in Nr uc (E) are of an irregular nature. These irregular fluctuations are exactly 
those in which one is interested in the theory of quantum chaos. 

To compare the quantal spectra of different systems, one has to get rid off the system- 
dependent constants in N(E), which is achieved by unfolding the spectra [23] by 

x n := N(E n ) , n = 1,2,3,... . (5) 

The unfolded spectrum {x n } has by construction a unit mean level spacing. In the sequel we will 
assume that the spectra are already unfolded, and in order to keep notation simple, N(x) will 
denote the spectral staircase of the unfolded energy spectrum, and Nr uc (x) the corresponding 
fluctuating part. Thus N(x) = x + Nr uc (x). 

To investigate the statistical distribution of the spectrum, several statistical measures have 
been introduced [23, 24]. The simplest quantity is 8 n , which is obtained from the mode fluc- 
tuations, i.e., the fluctuating part of the spectral staircase evaluated at the unfolded energy 
eigenvalues x n 

8 n := N (x n ) - N(x n ) = n - - - x n , (6) 

where N (x) := lim^o N ( x+e )+ N ( x ~ e ) . The quantity S n can serve as a good measure for the 
completeness of a given energy spectrum. 

Another statistic is the level-spacing distribution P(k, L) which is the probability density of 
the distances x n+ k+i—x n between energy levels. The well-known nearest-neighbour level-spacing 
is obtained for k = 0, P(L) = P(0, L). 

The number variance allows one to investigate medium- and long-range correlations of the 
spectrum. It is defined as the local variance of the number n(L } x) := N(x + L) — N(x) of 
unfolded energy-levels in the interval [x, x + L] 

E 2 (L,x) := ([n(L,x) - L] 2 ). , L>0. (7) 

For the averaging there are a priori different possibilities. We choose the following rectangular 
averaging 

cx 

</(*)>* = TTTTji / dx ' ( 8 ) 

where c > 1. The expectation value of n(L, x) is given by 

_^ cx 

(n(L,x))z = — / { N(x + L) - N(x) } dx 

(c — i)x J 

x 

^ cx 

= L + 7 / { N ^c(x + L)- N nuc (x) } dx , 

(c — l)x J 

x 

which tends to L in the semiclassical limit of x — > oo } since A^fl uc (a;) has zero mean. 

A related statistic, which also measures two-point correlations, is the spectral rigidity [25] 

A 3 (L,x) := ^min-j- j dy [N (x + y) - a - by] 2 ^j , (9) 

~~ X, 
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which is the average mean square deviation of the spectral staircase function from the best 
fitting straight line over an energy range of L mean level spacings. 

It has been conjectured that the energy level statistics of integrable systems can be described 
by a Poissonian random process [26], whereas the statistics of classically strongly chaotic systems 
should be reproduced by the corresponding RMT distributions of the GOE or GUE [4, 5, 6]. 
It turned out that this is only true for short-range correlations, e.g., measured by the nearest- 
neighbour level-spacing distribution P(0, L) [1, 27]. A convenient measure for the long-range 
correlations is provided by the number variance E 2 (L, x), which, however, does not discriminate 
the different short-range behaviours since one gets the behaviour E 2 (L,x) = L + 0(L 2 ) solely 
from the fact that the spectra are unfolded. The RMT results for large L are Eq OE (L) ~ 
^ln(2^L) and S^ UE (L) ~ ^ln(2^L) which is in strong contrast to the observed saturation 
of the number variance (see, e.g., [7, 11, 28, 29, 30]). The first chaotic system for which the 
saturation behaviour was studied was the Hadamard-Gutzwiller ensemble [31]. Berry has given 
semiclassical arguments in favour of the saturation behaviour of the number variance [30, 32] 
and of the spectral rigidity [2]. It turned out that with L max := ay/x~, where a > is a 
system-dependent constant, one has the universality regime for L <C iv max , in which the spectral 
statistics are expected to be described by random matrix theory. For L ^> L mSLX the number 
variance oscillates around its saturation plateau E 2 (x). This saturation plateau increases with 
increasing energy x in the limit of large x. 



Ill E(k, L) functions 



The E(k } L) functions, k = 0,1,2,..., are defined as the probabilities to find k levels in a 
randomly chosen interval of length L [25, 33]. In RMT the E(k, L) functions do not depend on 
the considered interval. For real systems, however, this stationarity is not fulfilled. Therefore, 
one considers E(k, L } x), the probability of n(L, x) being equal to k, where x is chosen uniformly 
distributed in the interval 

_^ CX 

E(k } L } x) = — — — / 4, n (L,^) dx , (10) 



(c — l)x J 

x 

where h n , m is the Kronecker symbol. Since E(k } L } x) and E 2 (L,x) are only directly related if 
one chooses the same averaging in both statistics, we consistently use the rectangular averaging 
(8). In order to keep the notation simple, we will write E(k, L) and E 2 (L). 

From the definition follows that the E(k, L) functions have their maximum approximately at 
L = k since only unfolded spectra are considered. Furthermore, they satisfy the normalization 

oo 

E E(k,L) = 1 • (11) 

k=0 

The number of levels n(L) in an interval of length L has the expectation value 

oo 

HD) £ = £ kE(k,L) ~ L . (12) 

k=0 

The last asymptotic equality holds strictly in the case of RMT. 
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The E(k, L) functions determine the other statistics as we would like to illustrate in some 
examples. Denoting by F(k, L) the probability to count k levels in an interval of length L, 
which starts at an arbitrary quantal level X{ } one has the relation 

F(k,L) = £ E(i,L) . (13) 

The level-spacing distribution P(k, L), which is the probability to obtain exactly k levels in an 
interval of length L beginning at an arbitrary quantal level X{ and ending at the quantal level 
x i+k+1 , is given by 

P%L) = £ F(i,L)= ^ j: (k-i + l)E(i,L) . (14) 

The nearest-neighbour level-spacing is the special case k = 0. As a last example we would like 
to recall that the number variance S 2 (L) is given by 

oo 

E 2 (L) = J2(k-L) 2 E(k,L) . (15) 

k=0 

In the case of a Poissonian random process one has an analytic expression for the E(k, L) 
functions 

E Poisson {k,L) = ^L k e~ L . (16) 

In figure 1 the Ep i sson (k, L) functions are shown in comparison with the E(k, L) functions 
obtained from the first 100 000 quantal levels of the circular billiard being an example of an 
integrable system (see section VII). For small values of L there is quite good agreement of the 
E(k, L) functions with Ep i sson (k, L), but for larger values of L clear deviations are visible which 
reflect the saturation of the E(k, L) functions. 

III.l RMT behaviour for the Ep(k,L) functions 

RMT starts with the joint probability measure Vp(x\, . . . , xjv) dx\ . . . dx^ [25, 34], which is the 
probability of finding a level in each of the intervals X{ + dxi], i = 1, . . . , N } regardless of 
the labelling, i. e., Vp(x\, . . . , xjv) is symmetric under the interchange of X{ with x r For the two 
commonly used ensembles, the Gaussian orthogonal ensemble (GOE) and the Gaussian unitary 
ensemble (GUE), RMT yields [25] 



Vpixt, . . . ,x N ) = Cpexp \--J2 x l( nixt-Xjf ( 17 ) 




with (3 = 1 for the GOE and (3 = 2 for the GUE, where Cp is a normalization constant. Since 
RMT deals with matrices of infinite dimension the limit N — > oo has to be considered. 

The E(k, L) functions can be obtained from Vp(x\, . . . , xjv) in the following way. To get 
exactly k levels in an interval / = [X } X + L] one has to integrate Vp(x\, . . . , xjv) for the first k 
levels over the interval / and the remaining levels over the range excluding / 

E fj (k,L) = — — / dxi...dx k I dx k+1 ...dx N V / 3(x 1 ,...,x N ) , (18) 

(1\ — k)\ Ji Jr\i 
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Figure 1: The E(k } L) functions are shown for the Poisson distribution and for the circular 
billiard discussed in section VII where the quantal levels up to the 100 000th are used. For a 
given k, the corresponding E(k, L) function is given by the curve having its maximum at L ~ k . 



where the combinatoric factor is due to the symmetry with respect to the interchange of X{ with 
Xj. This function depends on the choice of the interval /, i. e., on X. If, however, the spectrum 
is stationary, i.e., invariant under translations, as it is the case for GOE and GUE [35, 36], the 
result is independent of this choice after N — > oo. For the systems which will be discussed later 
an energy dependence is observed, and thus the results depend on the choice of /, i. e., on X. 

This integration process leads to the study of the continuous, real, symmetric integral kernels 
[25, 33, 37] 

sm(^(x-y)) 

K(x,y) := - and K±(x,y) := K(x, y) ± K(x, -y) , x,y£[-l,l] , 

tt(x - y) 

, (19) 

since the Ep(k } L) are related to the Fredholm determinants T>(L^f) = det(l — "fK) and 
Z>±(L,7) = det(l — "fK±). For the GUE case (/3 = 2) one obtains 



k\ d~f k 

Thus the Fredholm determinant T>(L^f) can be expressed as 



(20) 

7=1 



D(L, 7 ) = E^-Tl'^M) • (21) 

k=0 



6 



The GOE case (/3 = 1) is obtained by denning analogously to (20) functions E±(k } L) from 
which the E\(k, L) are determined by the relations 

E+(k,L) = E 1 (2k,L) + E 1 (2k — 1,L) , k > 

E-(k,L) = E 1 (2k,L) + E 1 (2k + 1,L) , k > (22) 

with £i(0, L) = E + (0, L). The Fredholm determinants are computed by diagonalization of the 
integral kernel K(x,y) leading to the problem of solving the eigenvalue problem 

j\ dy K(x, y) fi(y, L) = X t (L) fi(x, L) , (23) 

where the eigenfunctions fi(x, L) are given by the prolate spheroidal functions. They possess a 
definite parity fi(x, L) = ±/ 4 -( — x, L) because of the parity symmetry (x, y) — > (—x } —y) of the 
integral kernel. The Fredholm determinants are then given by the products 

oo 

W,7) = Ila-TA^)) (24) 

8 = 

and 

oo oo 

2>+(L, 7 ) = n(l-7A 28 (£)) , 2?_(L, 7 ) = n(l-7A 2 ,--iW) • (25) 

8=0 8=1 

With the above equations it is possible to compute the Ep(k } L) alone from the eigenvalues Xi(L) 
of the prolate spheroidal functions. It is worthwhile to note that the original version provided 
in [33] requires the computation of integrals over the prolate spheroidal functions fi(x } L) } a 
procedure which is now obsolete. The differentiation has to be carried out analytically to 
obtain numerically stable results. For the GUE case one gets 

oo 

E 2 (k,L) = J2 A n (L) ... K k {L) II (l-A.W) • (26) 

h<-<ik „ r 3= ° . , 

jg{'l>-,'fc} 

For the GOE case one has similar expressions for E + (k } L) and E_(k } L) where only even and 
odd eigenvalues occur, respectively. The E\(k, L) are then given by (22). Although the product 
runs over infinitely many eigenvalues, a fast convergence is numerically observed since almost 
all eigenvalues are very close to zero or to one. The above representation of the E(k, L) without 
the integrals over the spheroidal functions has also valuable applications as the determination 
of power series in L [37, 38] or the asymptotic behaviour L — > oo of the Ep(k } L) [39], which is 
also investigated by other methods in [40, 41]. 

Let us now turn to the integral equation (23) which will be solved separately for even and 
odd parity. The prolate spheroidal functions are expanded in Legendre polynomials with fixed 
parity up to order 2N or 2N + 1, i. e., for even parity fi(x, L) = Yl^o o-k{L) P2k{%), an d f° r °dd 
parity fi(x, L) = ^^_ (ik{L) P2k+i(x)- To discretisize the eigenvalue problem, N + 1 points xi 
in the interval [0, 1] are chosen which provide N + 1 equations by inserting these points into 
(23). Let a := (a^), then the generalized eigenvalue problem Aa = XBa is obtained whose 
matrix elements for even parity read 

Ajk = J dy K(xj,y) P 2 k(y) and B jk = P 2k (xj) . (27) 
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Figure 2: The E(k, L) functions are shown for the GOE (full curves) and for the GUE case 
(dotted curves). 

Analogous matrix elements yield the eigenvalue problem for odd parity. The numerical solution 
of the eigenvalue problem provides the eigenvalues Xi(L) as well as the prolate spheroidal func- 
tions, which are determined by a. We choose N = 20, . . . , 25 which allows sufficiently accurate 
results. In figure 2 the E(k, L) functions are displayed for the GOE and the GUE case. 

Often in the statistical analysis of energy spectra, main emphasis is laid on the investigation 
of short-range correlations by considering, e.g., the nearest-neighbour level-spacing distribution 
P(0, L). The "corresponding" E(k, L) function is £"(0, L) which is the probability of finding no 
eigenvalue in an interval of length P, and thus it is also called gap-probability. In this case good 
approximations for P(0, P) can be obtained using (c.f. equation (14)) 

P(0,P) = — P(0,P) (28) 

from the approximations to P(0,P): In the case of GOE the nearest-neighbour level-spacing 
distribution P(0, P) is very well described by the Wigner distribution 

P GOE (0,P)^Pexp{-^P 2 } , (29) 

whereas for GUE one has 

P GUE (0,P)^P 2 exp{--P 2 } . (30) 

Using (29) and (30) in equation (28), one calculates P(0, L) by integrating twice. The integration 
constants are determined by the fact that one has to fulfill P(0, 0) = 1 and limj^oo P(0, L) = 0. 
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One obtains 



E 1 (0,L) ~ erfc [^Lj (31) 
E 2 (0,L) ~ exp ( ——L 2 ^j — Lerfc (-y=Lj . (32) 



Remarkably, the approximations for £i(0, L) and £"2(0, L) are rather accurate (absolute error less 
than 0.005), despite the fact that they are based on approximations, which then were integrated 
twice. 



III. 2 Phenomenological description of the E(k, L) functions 

In the case of GOE and GUE [31] as well as for classically chaotic systems [31, 42] it is found 
that the E(k, L) functions are excellently approximated for some k > k by 

E(k,L) ~ \ T ^ ^(- { -T7Tv] ■ ( 33 ) 



The relation of the width a(L) to the number variance S 2 (L) [31], 

a(L) = T?{L) , L->k , (34) 

is obtained by inserting the approximation (33) for k > k in eq.(15) and by approximating 
the sum by the Euler-MacLaurin summation formula. A numerical comparison shows that 
the Ep(k } L) functions belonging to GOE and GUE are well described by (33) already for 
ko = &RMT — 6- Since for chaotic systems it is observed that the E(k, L) functions agree with 
the Ep(k } L) functions for k <^ k s ~ L mSLX with k s increasing with energy, one gets for chaotic 
systems k ~ 6, in general. 

The approximation (33) describes a Gaussian distribution with width a(L) if the E(k, L) 
functions are considered as a function of k with L fixed. However, this is a discrete statistic 
since k is an integer. On the other hand, if a(L) = S 2 (L) is only slowly varying in a region of 
order a(L) around L it is justified to replace a(L) by a{k) in (33) 

E(k,L) ~ -T^^(- iJ ^f] , (35) 

i.e., a Gaussian behaviour for E(k } L) considered as a function of L with k fixed, which is a 
continuous quantity. 

This replacement is justified for GOE and GUE where one has for L ^> 1 



E^ OE (L) ~ — 

7T 



TT 2 



ln(27rL) + 7 + 1 - — 

o 



1 



Kve(L) * — Pn(27rL) + 7 + 1] , 

TT 

with 7 = 0.5772 . . . being Euler's constant. Because of this logarithmic behaviour the Ep(k } L) 
functions can be considered as Gaussians in L. 
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In the case of classically chaotic systems the E(k } L) functions saturate [31, 42] for some 
k > k s with k s being a system and energy dependent constant proportional to L mSLX . Saturation 
in the E(k, L) functions means that their width does not further increase but rather oscillates 
around a plateau value. Furthermore, the E(k, L) functions can be well approximated by the 
Gaussians (35) for k > min(A; s , £;rmt)- A higher amplitude of the E(k, L) functions corresponds 
to a smaller width a{k). The saturation of the width a(k) around a plateau corresponds to 
the well-known saturation of the number variance 

«oo = , (36) 

which again follows from relation (34). 

It is worthwhile to note that in RMT S 2 (L) and A 3 (L) are related by 

A^(L) = jj (£ 3 - 2L 2 r + r 3 ) E^ MT (r) dr . (37) 
o 

Assuming this relation to hold also for an individual system, one can derive [31] that the sat- 
uration value of the number variance is related to the saturation value Aoo of the spectral 
rigidity by 

- SAo. . (38) 

Numerical results [11, 28, 31] for strongly chaotic systems suggest that this relation is indeed a 
good approximation. 

After having discussed the consequences of a Gaussian behaviour for the E(k, L) functions 
for the number variance E 2 (L), we would like to turn to the spectral form factor F(t), which 
in RMT is connected with the E(k, L) functions by the relation [29] 

F(t) = 2ttt / dL sm(2TTTL) — E 2 (L) . (39) 
Jo dL 

Inserting the approximation (35) in (39) one approximately obtains for L ^> k 0} i. e., for r — > 0, 

oo 

F(t) ~ -4^ 2 r 2 J2 a ( k ) i 1 -^ 2 T 2 a(k)) e^ 2 T * a ^ cos(2ir kr) , (40) 

k=0 

which reproduces the small r behaviour of F(t) in the case of Poisson and RMT spectra. 



IV The distribution of r]i(x) 

Before we discuss the distribution of the mode fluctuations W(x) } we turn to the distribution of 

,I(I):= 75^ ' (41) 

where x is randomly chosen from the interval [T, cT] with c > 1, and L > being a fixed 
parameter. The distribution of t]l(x) with x, L — > oo is in a certain limit related to W(x). 
Because of the expectation value (n(L } x)) ~ L and the definition of E 2 (L,x), the distribution 
of t]l(x) for T — y oo has zero mean and unit variance. In the case of a numerical analysis 
with a limited number of quantal levels the distribution of t]l(x) is not smooth and displays 
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pronounced equidistant peaks since for fixed L the quantity n(L, x) — L is an integer minus a 
fixed number, and E 2 (L, x) is only slowly varying. A smooth distribution can only be expected 
in the limit L,x — > oo because of the increase of E 2 (L,x). There are rigorous results for the 
distribution of t]l(x) in the case of RMT [43] and for some integrable systems [3], see below. 

Let us briefly recall some facts about the characteristic function, moments and cumulants of 
a distribution (see e.g. [44, 45]). The generating function F(t) of a distribution function P(x) 
is defined by 

/•oo 

F(t) := / t x dP(x) . 



For the special choice t 
obtained 



which exists, since 



u 



real, the characteristic function /( 



F 



dP{: 



(42) 

u ) of the distribution is 
(43) 



dP(s 



< 



dP{: 



rriQ 



The characteristic function f(u) determines uniquely the distribution P(x). By setting t := e p 
fi real, in eq.(42) the moment generating function M(/i) is obtained 



M(/i) 



F 



dP(s 



E 

k=0 



k\ 



(44) 



if the integral exists. Notice, that the convergence is not ensured in contrast to the characteristic 
function (43). Here m k := Jfj^ x k dP(x) are the moments of P(x). 

The moments of the distribution of n(L, x) can be determined from the E(k, L) functions 
by the moment generating function 



M(/i) 



E E(k,L)e 

k=0 



E 

k=0 



k\ 



(45) 



whenever the series are convergent. For a centered Gaussian with variance a 2 the odd moment 
vanish, and the even moments are given by 



(2a 



2\k 



m 2 k 



r(fc + i/2) 



(2*Q! 2k 
k\2 k 



(46) 



The even moments are increasing with increasing k. If one deals with a Gaussian distribution, it 
is thus more desirable to use the cumulants instead of the moments, since the cumulants c k have 
the advantage that for a Gaussian distribution only the first two cumulants are non-vanishing. 
The cumulants c k , k > 2, have the additional advantage that they are translational invariant 
under a shift in the variable L in contrast to the moments, i. e., the origin, with respect to which 
they are measured, plays no role, except for c\ = m\. Thus, the cumulants c k , k > 2, for n(L, x) 
and for n(L, x) — L are identical. If one transforms a distribution into the so-called standard 



measure, i. e., with zero mean and unit variance, one has 



d := and C k 



Ck 
k/2 



(47) 



The standard measure of n(L } x) is t]l(x), of course. To obtain the cumulants C k of t]l(x), we 
first compute the cumulants c k of n(L } x) } and then transform them according to (47). If the 
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cumulants Ck with k > 3 vanish for L — > oo, the distribution of t]l(x) is a Gaussian in the weak 
limit. The cumulants have, however, the disadvantage that they cannot be computed directly 
from the data, i.e., with integrative or summatory processes. They can be computed from the 
moments or from the cumulant generating function, respectively, which is related to the moment 
generating function. The cumulant generating function reads 

CV) = In MM = £ c k ^ , (48) 

k=i K - 



so that the cumulants are given by 



(49) 



In RMT the asymptotic behaviour for L — > oo of the cumulants Ck for the distribution of 
n(L } x) is studied in [43]. Because of the stationarity of the random matrix ensembles one 
sets the energy variable x = — L/2 and considers the probabilities with respect to the interval 
[— L/2, L/2]. The second cumulant is c 2 (L) = In L + 0(1) for the GUE case and, therefore, 
the &th cumulant of t]l will be Ck — Ck/(-^i lnL) k / 2 . Thus a Gaussian behaviour is displayed if 
Ck — > for k > 3, i.e., Ck = o ((-^ In L) k ^ 2 ^j for k > 3. This behaviour is shown to be true for 
GUE, and a similar behaviour is obtained for GOE and GSE in [43]. Therefore the distribution 
of r]L is Gaussian in at least these three RMT ensembles in the weak limit for L — > oo. 

For integrable systems with eigenvalues E nufl2 = I[n\ — «i,n 2 — a 2 ) } where I(xi 7 x 2 ) is 
a smooth homogeneous function of second degree and («i,a 2 ) is a point in the unit sphere, 
the distribution of n(L } x) is studied in [3, 46]. It is proven that N(x + L) and N(x) are 
asymptotically independent for L/x — > and Lj^fx — > oo. This limit can be realized by the 
choice x — > oo with L = zx s , | < S < 1. 

Furthermore t]l(x) converges to a difference of two independent identically distributed ran- 
dom variables, whose limit distributions are the same as for W(x). In the limit x — > oo with 
L/y/x — y z > it is proven that t]l(x) has a non-Gaussian limit distribution. However, for 
z — y this distribution tends to a normal distribution. In addition rigorous results concerning 
the behaviour and saturation of the number variance can be found in [3, 46, 47]. 

For the chaotic quantum systems discussed later, a Gaussian behaviour is observed for the 
E(k } L) functions for sufficiently large k values as a function of L, see eq.(35). Due to the 
saturation this implies that the approximation (33) is fulfilled for large L. Using this approx- 
imation in (45) and (48) one obtains a series, which is approximated by an integral using the 
Euler-MacLaurin summation formula. Evaluating the integral yields for large L 

C(fi) ~ Lfi + ^^/i 2 , L»l . (50) 

Therefore one obtains c\ — > L } c 2 — >■ a(L) and Ck — >■ for k > 3 which yields Ck = 0, k > 3, of 
course. Thus the distribution of t]l = (n(L) — L) / yj 'a(L) is a normalized Gaussian in the weak 
limit L — y oo ; if the Gaussian behaviour is established for E(k,L). 

In the case of a Poissonian random process one obtains with formula (16) for Ep i sson (k } L), 
using (45) and (48), the cumulant generating function 

C(fi) = L(e^-l) . (51) 

Together with eq.(49) it immediately follows that for all cumulants Ck = L in the case of a 
Poissonian random process. Since C'k = L/L k / 2 — > 0, k > 3, for L — Y oo, the distribution of t]l 
approaches a Gaussian for a Poissonian random process in the sense of a weak limit. 
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V Mode-fluctuation distribution P(W) 



As already discussed, the description of spectral statistics using finitely many energy levels 
by means of RMT is restricted to short-range correlations because of, e.g., the saturation 
of the number variance and rigidity, contrasting the expected logarithmic increase in RMT. 
Furthermore, even for short-range correlations deviations were found for a class of strongly 
chaotic systems, the so-called arithmetical systems [9, 10, 11, 13]. Therefore in [14, 15] a 
different statistical measure as signature of quantum chaos is proposed which is given by the 
probability distribution of the suitably normalized fluctuations of Nr uc (x). 

20 I . . 1 1 1 



15 - 




-15 - 



-20 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

20000 40000 60000 SS0000 x 100000 



Figure 3: Nr uc (x) for the integrable circular billiard. 

In order to define this statistical quantity, some considerations in advance are necessary. 
Nft ac (x) is shown in fig. 3 for the integrable circular billiard (see sec. VII. 2), using the first 
100 000 eigenvalues of the even symmetry class with respect to reflection at the diameter. Since 
Vfl U c(s) behaves very irregularly one might consider Nr uc (x) as a random function of x and look 
for its statistical properties. The increase in the amplitude of the fluctuations is clearly visible 
in fig. 3. Therefore, one defines the variance 

D ^ := ( c -l)T } [Nn ^ ]2 dy ' < 52 ) 

with c > 1, where S(c) is a correction to be explained below. With this one can construct the 
following random variable 

W(x) := %M , (53) 
'D(x) 
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where x is randomly chosen from the interval [T, cT]. This statistic was first investigated in the 
special case of the eigenvalues of the Laplacian on a torus [27]. A similar statistic for the GUE 
without the normalization was studied in [48]. If a limit distribution for T — > oo of this random 
variable W(x) exists, it has by construction a second moment of 1 if the second moment exists. 
Furthermore, the first moment of this limit distribution is zero because by definition Nr uc (x) 
describes the fluctuations of N(x) around the mean behaviour N(x) = x. 
The conjecture proposed in [15] can now be formulated as follows: 

For bound conservative and scaling systems the quantity W(x), eq. (53), possesses a limit 
distribution for T — > oo. This distribution is absolutely continuous with respect to the Lebesgue 
measure on the real line, with a density PiW) 



cT 



t 1 ™ (c-l)T J giW{x)) P{X/T) dX= I giW) PiW) dW ' (M) 

^ X -oo 

where g(x) is a bounded continuous function, and pit) > is a continuous density on [l,c] with 
^ 1 f 1 c p(t)dt = l. 

Furthermore, the limit distribution has zero mean and unit variance, 

oo oo 

j dWW P{W) = , j dWW 2 P{W) = 1 . (55) 

— oo — oo 

// the corresponding classical system is strongly chaotic, having only isolated and unstable peri- 
odic orbits, then P(W) is universally a Gaussian, P(W) = y= exp( — ^W 2 ) . In contrast, a 
classically integrable system leads to a non-Gaussian density PiW). 

In the case of several integrable systems eq. (54) is proven, and furthermore, it is shown 
that the corresponding limit distribution decays faster than a Gaussian. These results will be 
summarized in the following section. 

Some comments are in order: 



First notice that eqs. (55) do not follow from (54) by choosing g(x) = x and g(x) 



x 



2 



respectively, since g is supposed to be bounded. The question of the existence of higher moments 
is rather non-trivial, as will become clear in the case of integrable systems in the following section. 

Secondly, if a chaotic system has non-isolated periodic orbits, the conjecture can be general- 
ized by taking their contribution N\ ong (E) into account in "Weyl's law" N(E), eq.(4). Unfolding 
the energy levels in this way yields a Gaussian behaviour for the distribution of W(x). Notice 
that in the case of other long-range contributions which vanish for x — > oo, it is sometimes 
necessary to take them into account in N\ ong (x) in order to approximate the limit distribution 
sufficiently well with only a finite number of energy levels. 

Thirdly, note that the limit distribution PiW) is independent of the choice of the probability 
density p(t). This is easily established, if one can prove eq. (54) in the case of pit) = 1, by 
approximating an arbitrary p(t) £ L 1 ([l,c]) by step-wise functions (see [49, 50]). By choosing 
different densities pit) the distribution of x on [T, cT] is p(x /T) dx /T , which enables one to 
make a change of variables [3, 47]. E.g., one can choose pit) = 1, which corresponds to an 
equidistribution of x on [T, cT]. Choosing pit) = 2 (^/Z-i) ^~ 1 ^ 2 corresponds to averaging with 

respect to d(x 1 ^ 2 ), i.e., yjx is distributed uniformly on [y/T, a/cT]- Similarly for p(t) = -^t 
one obtains an equidistribution of x 2 on [T 2 , (cT) 2 ]. 

Berry's semiclassical analysis [2] provides predictions of the asymptotic behaviour x — > oo of 
the spectral rigidity A 3 (L } x) } i.e., about the energy dependence of the saturation value A O0 (x) 
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for x — y oo. Analogous arguments carry over to the energy dependence of the variance D(x). 
For generic integrable billiards one obtains 



D(x) ~ Clnty/X , (56) 

where c\ nt is some non-universal constant which depends on the considered system as well as on 
the chosen average. In the case of generic classically chaotic systems with anti-unitary symmetry 
(e. g. time-reversal symmetry), these arguments yield 

D{x) ~ J-lnz + C , (57) 

Z7T 



and for systems without such a symmetry 



D(x) ~ -^-Inx + C . (58) 
4tt 2 



In the case of arithmetical chaotic systems arguments based on the diagonal approximation of 
the spectral form factor lead to [51] 

D(-) - , (59) 

where A is the area of the corresponding arithmetic surface. However, the arithmetic hyperbolic 
triangle discussed in section VII. 1 possesses a behaviour which is more conform to c\y/x + c 2 
like integrable systems [11]. These results can be used in asymptotic considerations as analytic 
expressions for the variance D(x) (see section VII. 3 for details concerning numerical tests of 
the conjecture). 

We are now in a position to explain the factor S(c) occurring in the variance D(x) in eq.(52). 
To obtain a distribution of W(x) with variance one it is necessary to normalize Nr uc (x) with 
the variance at x. However, with D(x) := ([Nr uc (x')] 2 ) x one uses a variance corresponding to 
a point x 1 £ [x,cx] for the normalization which differs from the local variance at x for large 
averaging parameter c or for sufficiently increasing D(x). Therefore we now study the alteration 
caused by the averaging (8) in the functional dependence on D(x). 

Considering a power law for D(x) = ax 1 which includes (56) with 7 = f, one obtains 

D(x) = (D(x')) £ = D(x) (60) 

which in turn leads to 

H W - J^I) fM -' = 5 ' (61) 

For a functional dependence D(x) = alnx + 6, covering the chaotic cases (57) and (58), one 
obtains 

D(x) = D(x) + a ^ - 1 ^ . (62) 

Thus for chaotic systems (excluding arithmetical chaos) one gets S(c) = 1. Asymptotically the 
constant plays no role, but for numerical computations it is sometimes necessary to take this 
additive constant into account. 

The above formulation of the conjecture is slightly different from the one given in [15]. 
Therefore let us now describe in short the connection between both formulations: In [15] the 
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mode-fluctuating distribution was introduced for bound conservative scaling systems with at 
least two degrees of freedom. There the variable x = E a , with a being the scaling parameter, was 
used as the basic variable. From a theoretical point of view this choice seems more appropriate 
because for scaling systems the natural variable is x in Gutzwiller's trace formula [7]. For 
numerical considerations, however, it is much better to have an equidistribution on the whole 
unfolded energy interval, which gives every level the same "weight". For example, a uniformly 
distributed random variable in p pronounces the contribution of the lower energy levels, which 
is not desired. Nevertheless, there is no contradiction between both formulations, since both 
are equivalent, which can be shown by appropriate choices of the density pit): 

Choosing p(t) = "^"^ t"" 1 , one can change from the x formulation to the formulation in the 
energy variable E. The reverse direction is obtained by p(t) = jjj^zTj^t^ -1 . If one neglects for 
large energies E the circumference and the constant term in the generalized Weyl's law (4) it is 
trivial that the formulation in E and the one in the unfolded variable x are equivalent. 

For classically chaotic systems a semiclassical approximation for Nr uc (x) is given by means 
of Gutzwiller's periodic-orbit sum. Here Nr uc (x) can be expressed using the dynamical zeta 
function Z(s) which involves only properties of periodic orbits like their lengths, stabilities and 
Maslov indices. The non-trivial zeros of the dynamical zeta function Z(s) are directly connected 
with the quantal levels of the corresponding system. Using the argument principle, one obtains 

^fluc(p) = -arg Z(ip) . (63) 

TT 

The tilde is chosen to distinguish the different arguments, i.e., the momentum p instead of the 
unfolded energy x. 

An interesting example is provided by the non-trivial zeros of the Riemann zeta function C( 5 )- 
Their statistics behave in many respects like the eigenvalues of some hypothetical classically 
chaotic system without anti-unitary symmetry. Here one has 

iVfluc(p) = -arg C (7: + ^ 

TT \Z 

It can be shown that for a certain class of ^-functions, which includes the Riemann zeta function, 
the corresponding W(p) = N ^ c ^ has a Gaussian limit distribution [52, 53], see also [54]. 

^/lnlnp 

In order to shed some more light on the meaning of the conjecture, we emphasize the fact that 
a Gaussian limit distribution reflects the validity of a central limit theorem for the fluctuations. If 
one considers all probability distributions with density and with unit variance, then the spectral 
entropy 

+ oo 

£[P] ■=- f dW P{W)\nP{W) 

— CO 

is maximized for a Gaussian distribution with mean zero. Since £[P] measures the mean 
unlikelihood for W(x) to be of a specific value, it provides a quantitative measure for spectral 
randomness. Thus the content of the conjecture is that classically strongly chaotic systems 
have maximally random quantum spectra. We should mention that from the numerical point 
of view £[P] is not a suitable quantity to distinguish between a Gaussian and a different limit 
distribution [42]. 

Before we present the relation between t]l(x) and W(x) } we need a relation between the 
number variance E 2 (L,x) and D(x). To be consistent with the definition of the variance D(x), 
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eq.(52), we use an analogous definition for E 2 (L,x), which gives the number variance at x and 
not at some x £ [£, cx] as explained above. Such a definition differs from (7) only for systems 
with a non-logarithmic behaviour of D(x). Thus we get 



E 2 (L, x) = , A / dx [N(x + L)- N(x) - Lf 
(c — l)x J 



dx [N auc (x + L) - N auc (x)] 



2 



(c — l)x J 

x 

= J dx {[NancXx + L)f - 2Nn uc {x + L)iV fluc (a;) + [iV fluc (a;)] : 

X 

If one considers the limit of large x and L in such way that one is in the saturation regime, 
i.e., L L max = Cy/x~, one expects A^fl uc (a; + L) and A^fl uc (a;) to become independent, which 
is proven in [3] for a class of integrable systems. Under this assumption the integral over 
^Vfluc(^ + -^)^Vfluc(^) vanishes. The last summand yields D(x), and one can show that the first 
summand yields D(x), too. To that end let us fix L = Cx s , with | < S < 1, and consider the 
first term for x — > oo 



dx [N nuc {x + L)f J dx [N nuc {x)} 2 -.jdx [N nuc {x)} 2 = ^— ^-D(x) . 

x £(1+C£«- 1 ) £ " 

For large x the upper bound tends to cx and the lower bound to x, because of 8 < 1. Thus we 
finally have 

E 2 (L, x) ~ 2 D(x) for large x and L = Cx\ \ <8 <l . (64) 

In comparison with t]l(x), the W(x) statistic has the advantage that one need not care about 
the different cases of x, L — > oo, because one is always in the non-universality regime L ^> L mSLX . 
Both statistics are related (the basic idea is due to Bolte [55]) in the following way: If one 
considers in 

. , N(x + L) — N(x) — L [N{x + L)-{x + L)]-[N{x)-x] 
r] L (x) = = (65) 

iyhTx) Jx*(l,x) 



the limit of large x in such a way that one is in the saturation regime by setting as above 
L = Cx s , | < 8 < 1, we can use (64) 



rj L (x) 



V2 



[N(x + L) - (x + L)\ [N(x) - x] 
D(x) \/D(x) 



(66) 



The second fraction in the bracket is W(x). If in the limit of large x, L with L = Cx\^<8<1 } 
the random variables N(x + L) and N(x) are statistically independent, the right hand side is 
the sum of two independent random variables. Assuming that t]l(x) is Gaussian distributed for 
L — y oo, then each of the random variables on the r.h.s. has a Gaussian distribution [44]. 

Furthermore, if W(x) has a Gaussian distribution, also the first term in eq.(66) has a Gaus- 
sian distribution. If both random variables are independent, t]l(x) is also Gaussian distributed. 
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The variance is one, because the variance of the sum of two independent random variables is 
the sum of each variance [44]. 

Thus for x,L — > oo with L = Cx s , | < S < 1, a Gaussian behaviour of i]l(x) implies a 
Gaussian behaviour ofW(x) and vice versa. 

Assuming the conjecture, it follows that in the above limit the energy levels of integrable 
systems cannot behave like the outcomes of a Poissonian random process, since a Poissonian 
distribution for the E(k, L) functions implies a Gaussian distribution of t]l, which has a Gaussian 
distribution for W(x) as a consequence. 

Furthermore, it should be remarked that the above discussion shows that W(x) is a statistic 
which measures fluctuations in the saturation regime, i.e., global fluctuations. 



VI Exact results for integrable systems 

Recently a lot of results concerning the limit distribution PiW) of W(x) in the case of integrable 
systems were obtained (summaries can be found in [46, 56, 57]). For several integrable systems 
(see below) it is possible to show that N(E) can be decomposed as follows 

N(E) = —E + E 1/4 Q(Ve) ■ 

47T 

The remainder term E l l 4 ®{\f~E) fluctuates around zero with an amplitude increasing as E l l 4 
and corresponds to Nr uc (E). The first term is Weyl's law. Since the considered systems have 
no boundary, there is no additional term involving C, and also the constant C is zero. 

The basic idea for proving the existence of a limit distribution of ®{\f~E) is to express N(E) 
in terms of a lattice-point problem and then to exploit the Poisson summation formula to obtain 
an expression of the remainder term as a trigonometric series. If the remainder term is an 
almost periodic function [58] (see also [49, 56]) then it has a limit distribution [49, 59]. To be 
more precise, we have to give some definitions (following the exposition in appendix II of [56]): 

The space B p of pth power summable Besicovitch almost periodic functions is defined as 



follows: / G B p , if for any e > there exists a finite set of numbers (a^u^), 
Uk G R, a k G C, such that 



k = 1, 



N(e) 



lim sup 

T^oo 



N(e) 

f(t) - a k exp(iuj k t) 



k=l 



dt 



i/p 



< e 



This means that every / G B p can be approximated in the mean arbitrarily well by trigonometric 
series. In B p one has the semi-norm 



bp 



lim sup 

T^oo 



y / \f(*)\ P 



dx 



i/p 



such that ||/ — g\\Bv = for /, g G B p defines an equivalence relation. For 1 < pi < p 2 one has 
B P2 C B P1 and for / G B P2 the inequality ||/||bw < holds. 

A formal trigonometric series J2T=i a fc ex p( — i^kt) is called generalized Fourier series of / 
with respect to B p ', if limjv->oo | \f(t) ~ SfcLi a k exp(— iu k t) \ \bp = 0. For a given / G B p its 
generalized Fourier series is well defined and unique. 
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Let us further define the space H } with / £ H } if for any e > 0, there exists a finite set of 
numbers (a^, u>k), k = 1, . . . , N(e) } u>k £ R, £ C, such that 



lim sup — / min | 1, 

X^oo T J 

o 



N(e) 

f(t) - a k exp(iuj k t) 
k=i 



dt < e . 



This assumption is weaker than / £ B 1 , since B 1 C H. 

In [49] it is proven using ergodic-theoretic and probabilistic methods that the condition / £ H 
is sufficient to establish that fit) possesses a limit distribution 



+ oo 



lim — 

X^oo T 



g(f(t)) dt 



g(x)u(dx) 



for every bounded continuous function g(x) on R. If in addition / £ B p } p > 1, then 



+ oo 



\x\ p u(dx) 



1/p 



< oo 



1 



+ 00 



-oo 



Thus for / £ B p the absolute moments up to the pth order exist. 

In a series of papers it is proven for several integrable systems that the corresponding nor- 
malized remainder term ®(\/E) belongs to H, or even B p . In these proofs N(E) is expressed in 
terms of a lattice-point problem. A very old problem of this type, due to Gaufi, is the classical 
circle problem, where one considers the difference between the number of lattice points inside a 
circle of radius R and the corresponding area irR 2 . For a long time the aim was to find uniform 
bounds of the remainder term. In [59] the study of limit distributions of the remainder term 
by means of almost periodic functions was started. From our point of view one is led to this 
problem if one considers the eigenvalues of the Laplace operator on a two-dimensional torus. 
These read (in suitable units) E n , m = n 2 + m 2 , n, m £ Z. Thus N(E) is given by the number of 
lattice points inside a circle of radius R = \[~E. Applying the Poisson summation formula gives 



'R 2 



N(R 2 



ttR 2 



E 

^ezl\{(o,o)} 



Ji(27r|n|i2) 



R 



where |n| = \jn\ + n\. Using the asymptotic behaviour of the Bessel function one obtains for 
the normalized mode fluctuations 0(i?) = Afl uc (i? 2 ) / ' \[~R 

, s 1 ^ cos(2dn|i?- , 1x 

n ez 2 \{(o,o)} l"l 

The normalized remainder term can then be rewritten as a sum over square free numbers m 

®{R)= E f m (V^R) + 0(R- 1 ) , 



square free TTi 



wi 



ith 



fm(t) 



3/4 ^ an cos ~ t 



71 E 

k=l 
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where gk, m '■= i^{ n £ Z 2 | \n\ = ky / m}. The above representation of ®(R) can be interpreted as 
a flow on an infinite-dimensional torus with incommensurate frequencies, where R is considered 
as the time-parameter. By this the problem of finding a limit distribution becomes a problem in 
ergodic theory. The above representation can also be interpreted as a trace formula, where the 
sum runs over the lengths \jn\ + n\ of primitive periodic orbits of the geodesic motion on a two- 
dimensional torus, and the sum in the representation of f m (t) corresponds to multiple traversals. 
It is proven in [59] that the corresponding normalized remainder term ®(\/E) belongs to B 2 
and possesses a limit distribution. Furthermore, the limit distribution is absolutely continuous 
with respect to the Lebesgue measure, i. e., it can be written as u(dx) = P(x) dx by means of a 
density P(x). In addition, estimates for derivatives of the density are given and the convergence 
of the (absolute) moments is proven up to the ninth moment. In [60] it is proven that the third 
moment is negative, which implies that the limit distribution is non-Gaussian. 

In [50] the more general case of the number of lattice points inside a shifted circle is con- 
sidered. It is proven that the corresponding normalized remainder term has an absolutely con- 
tinuous limit distribution and that the first and the second moments exist. The corresponding 
probability density P(x) decays faster than a Gaussian: for every e > and \x\ > x (e) 

< P(x) < exp(-|x| 4 - e ) , (67) 

which again shows that the limit distribution cannot be Gaussian. In addition in [61] it is proven 
that the variance of the remainder term is a wild function of the shift of the circle. 

In [49] the existence of a limit distribution and the first two moments is shown for the 
remainder term ®(\/E) of the number of lattice points inside a smooth strictly convex shifted 
oval including the origin. 

In [62] the existence of a limit distribution for the normalized remainder term of the shifted 
Laplace operator (zV — a) 2 on the two-dimensional torus T 2 (27rai, 27ra 2 ) is proven. It is shown 
for transcendental ratios ai/a 2 that the limit distribution possesses a density, which fulfills a 
slightly stronger condition than (67). In addition these results are extended to a general class 
of convex ovals. From these results it easily follows that for the integrable rectangular billiard 
with sides a and b and for the right-angled isosceles triangular billiard the corresponding limit 
distribution is non-Gaussian. 

In [63] the existence of the limit distribution is proven for the eigenvalues of the Laplace 
operator on certain surfaces of revolution. In particular, in this case N(E) can be written in 
terms of a lattice-point problem for the Einstein-Brillouin-Keller quantization. It is proven that 
the normalized remainder term is of Besicovitch class B 2 and its generalized Fourier series is 
given by a sum over closed geodesies. Under an additional assumption it could be shown that 
the limit distribution has a density, which fulfills < P(x) < C exp( — Ax 4 ), C, A > 0. Also in 
[63] a more general class of surfaces of revolution is studied, where the remainder term can be 
decomposed in the following way 

Nfi uc (E) = N(E) -4~ E = El/3 H^E) + E 1/4 F(VE) . (68) 

Here $(^E) is a finite sum of bounded periodic functions, which implies that its distribution has 
compact support. F{\f~E) is an almost periodic function of Besicovitch class B 2 , and thus has a 
limit distribution. In the above case the leading behaviour of Nr uc (E) is determined by the first 
term E 1 / 3 <&(\/E). Considering the behaviour of F{\f~E) corresponds to study the distribution 
of [Nq uc (E) — E 1 ! 3 <&(\/E)]/ E 1 ! 4 . This way of splitting is similar to the one introduced in eq. (4) 
involving the term N\ ons (E). 
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In [56] the remainder term for the eigenvalues of the Laplace-Beltrami operator on a two- 
dimensional torus with a Liouville metric is investigated. The geodesic motion on such Liouville 
surfaces can be seen as the most general case which is completely integrable. Using a refined 
Einstein-Brillouin-Keller quantization the study of the remainder term is reduced to a lattice- 
point problem. It is proven under a certain condition of non-degeneracy that the normalized 
remainder term is of Besicovitch class B l . For generic surfaces the distribution of ®(\/E) has 
a density, which decays faster than a Gaussian. In [64] a refined estimate for the decay of the 
distribution is proven. Also an explicit expression for the Fourier transform of the density by 
means of the Fourier amplitudes of ®(y/E) is given. In addition a trace formula for Liouville 
surfaces is derived. 

In [65] the distribution of eigenvalues of Schrodinger operators on manifolds all of whose 
geodesies are closed are studied. It is proven that the corresponding ®(x) is in B 2 } and that the 
limit distribution is an equidistribution on a finite interval. As a remark one should add that in 
the case of the one-dimensional harmonic oscillator and of a particle in a one-dimensional box 
the corresponding limit distribution is also a box-shaped function. 



VII Discussion of the considered systems 

VII. 1 Models for quantum chaos on surfaces of negative curvature 

The motion of a point particle on surfaces of negative curvature is of special historical interest 
in the context of chaos since already in 1898 with them Hadamard provides the first example 
of strongly chaotic systems [66]. In 1939 Hopf proves that negative curvature necessarily leads 
to ergodicity [67], and for the special case of constant negative curvature it is later proven that 
it leads to strongly mixing systems and its time-evolution operator has Lebesgue spectrum [68]. 
For compact Riemann surfaces of constant negative curvature Sinai proves that the systems 
belong to the class of A"-systems [69], i.e., they show hard chaos. However, the special role 
these surfaces play is not only of historical nature, since, in addition to the proven property of 
being a A"-system, they also possess due to the Selberg trace formula [70] an exact periodic-orbit 
trace formula. The latter is valid only semiclassically for general systems [7]. 

A convenient model for surfaces of constant negative curvature is the Poincare disc T> } 
consisting of the interior of the unit circle in the complex z-plane (z = X\ + ix 2 ) } endowed with 
the metric 

4 

9a = T\ 2 2v 5 a ' *>i = 1 > 2 • ( 69 ) 

(1 — x\ — X2) 

This metric has constant negative Gaussian curvature K = —1. This choice fixes the length 
scale. The classical dynamics is determined by the Hamiltonian 

1 ■ - dx^ 

To obtain a compact Riemann surface from the infinitely extended hyperbolic plane one tes- 
sellates the plane into copies of one fundamental domain by the action of a Fuchsian group. 
This can be viewed as cutting out a polygon from the plane in such a way that opposite edges 
can be glued together, i.e., periodic boundary conditions are imposed [7, 71, 72]. The simplest 
compact case allowed is a tiling polygon with eight edges, yielding a surface of genus two, which 
is topologically a sphere with two handles. For more details about these hyperbolic octagons 
see [7, 72]. 
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The set of all octagons provides an ensemble of /<f -systems, which can be investigated in 
the sense of RMT, i. e., an ensemble of systems with the same dynamical properties. The usual 
statistical analysis of a quantum chaotic system concerns in general only an individual system. 
In the case of the octagon ensemble, however, there is the possibility to study the statistical 
properties of a physical ensemble, where individual properties are averaged out. A statistical 
analysis of the E(k, L) distributions, and other important statistics of the octagon ensemble, is 
carried out in [31]. To our knowledge this is the only case where a statistical analysis of physical 
systems is carried out in the spirit of RMT. But since the main conjecture of quantum chaos is 
that the statistical behaviour of RMT is valid for a typical individual quantum chaotic system, at 
least at short-range correlations, we would like to study the properties of an individual octagon 
instead of a subset of the whole ensemble. 

The so-called regular octagon, which is the octagon with the highest possible symmetry, 
where all corners possess an angle of 7r/4, was the first being investigated more closely with 
respect to quantum chaos [8, 72, 73]. Because of the high symmetry one has to desymmetrize the 
regular octagon in order to compare its statistical properties with RMT. The one-dimensional 
irreducible representations of the symmetry group corresponds to a hyperbolic triangle with 
angles a = -|, (3 = -| and 7 = § having 1/96 the area of the regular octagon. There are 
eight possible combinations of Dirichlet and Neumann boundary conditions along the edges a, b 
and c lying opposite the angles a, (3 and 7, but only four provide subspectra of the regular 
octagon. One choice providing a subspectrum of the regular octagon is realized by imposing 
Dirichlet boundary conditions on all three edges, as considered in [72, 74]. Here we consider 
two particular combinations of Dirichlet and Neumann boundary conditions, of which one is 
contained in the spectrum of the regular octagon, called triangular billiard B in the following, 
and the other, called triangular billiard A, is not. Both have Neumann boundary conditions 
on side c, but billiard B has Dirichlet b.c. on side b and Neumann b.c. on side a, whereas 
billiard A has Dirichlet b.c. on side a and Neumann b.c. on side b. The statistical analysis of 
both billiards is discussed in [11]. The reflection group which generates triangle B belonging 
to the regular octagon is arithmetic and thus implies exponentially increasing multiplicities in 
the length spectrum of periodic orbits. Arithmetical quantum chaotic systems are now known 
to violate the RMT behaviour even for short-range correlations [9, 10, 11, 12, 13], and the 
triangular billiard B shows all features of arithmetic systems. The triangular billiard A and the 
other triangular billiards with boundary conditions which do not belong to the regular octagon, 
show no peculiarities in the quantal spectra although they have nevertheless the same degenerate 
spectrum of periodic orbits however with different characters attached to the periodic orbits. 
The following statistics are based on the first 1000 quantal levels for each triangle. 

In addition to these two triangles we also investigate an asymmetric hyperbolic octagon 
which only possesses one involutary symmetry which is always present in hyperbolic octagons. 
It is not arithmetical and can easily be desymmetrized with respect to the involutary sym- 
metry. The statistical analysis concerns the subspectrum belonging to the positive parity class 
comprising the first 1127 quantal levels. 

VII. 2 Euclidean billiards 

In this section examples of two-dimensional Euclidean billiards will be introduced. They are 
given by the free motion of a point particle inside a compact domain 0, with elastic reflections 
at the boundary. 
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cardioid billiard 



stadium billiard 




Figure 4: Cardioid billiard, stadium billiard, truncated hyperbola billiard and one of the iso- 
spectral billiards. 

1. Cardioid billiard 

The boundary curve of the cardioid billiard (see fig. 4) is given in polar coordinates by 
p(4>) = f+cos (f). The cardioid billiard is the limiting case of a family of billiards introduced 
in [75, 76], whose boundary is defined by a conformal mapping of the unit circle. For the 
cardioid billiard it is proven that it is ergodic, mixing, a /\-system and a Bernoulli system 
[77, 78, 79, 80]. Due to the reflection symmetry of the cardioid billiard, one can divide 
the quantum mechanical wave functions into odd and even eigenfunctions. In [28] no 
observable differences between the statistics of the odd and even energy spectra are found 
(this is also observed in [81]). Thus we use in the following analysis the hrst 10 000 energy 
levels for odd symmetry (for a detailed investigation of the spectral statistics for both 
symmetry classes see [28, 42]). The set of eigenvalues was calculated [82] by means of the 
conformal mapping technique [76, 83]. 

N(E), eq. (4), holds with A = C = 6 and C = for odd symmetry. There are no 
prominent slowly varying long range correlations observed [28], such that N\ ong (E) = 0. 

2. Stadium billiard 

The stadium billiard [18, 19] is given by two semi-circles of radius I joined with two 
parallel lines of length 2a (see fig. 4). In our calculations we choose a = 1 .8. It is proven 
that the stadium billiard is ergodic, mixing, a /\-system, and a Bernoulli system [18, 19]. 

We will consider the completely desymmetrized stadium billiard with everywhere Dirichlet 
boundary conditions, for which we calculated the hrst 6000 eigenvalues using the boundary 
element method. The constants in eq. (4) for a quarter stadium with Dirichlet boundary 
conditions are given by A = a+j, C = 2+|-+2a and C = The additional term N\ ong (E) 
is not zero because of the so-called bouncing-ball orbits. These are a family of neutral 
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periodic orbits bouncing perpendicularly between the straight lines. Their contribution 
to the spectral staircase is given by [20] 

*W£) = a ~ t ^ - ^ e (VI - v^F) - {^e - ±Ve) , 

where = ir 2 n 2 are the eigenvalues of a particle in a one-dimensional box of length 1, 
and is the Heaviside step function. 

There are further contributions from the bouncing-ball orbits due to edge effects [20, 21]. 
The integrated results of the formulas derived in [20] are 

— cos(2/E) 

1 ~ sin(2n/E) _ J ; \f~E~ = tt 

n " \ i[^modl]-| ; ^ tt " 

Both contributions are bounded and thus play no role for the limit distribution for E — > oo. 
But since we have only a finite number of energy levels, we find a slight improvement if 
we take them also into account in N\ ong (E). 

3. Truncated hyperbola billiard 

The truncated hyperbola billiard (see fig. 4) is a compact model of the unbounded hy- 
perbola billiard, which is studied in detail in [84]. The boundary is given by the positive 
x-axis, the diagonal y = x and the hyperbola y = 1/x, which is truncated at A = (a, I /a) 
by a circular arc in such a way that the first derivative is continuous and the arc is perpen- 
dicular to the x-axis. Using the boundary element method the first 1851 eigenvalues with 
E < 10 000 are calculated for the truncation parameter a = 7.1175806, which corresponds 
to the value chosen in a recent microwave experiment [85] . 

The spectral staircase function displays a prominent long-range modulation, which can 
be explained by an additional term in Gutzwiller's trace formula [22]. This contribution 
arises from families of closed, non-periodic orbits which are reflected in the point A, 
where the curvature is discontinuous. This contribution tends to zero for E — > oo, but 
is necessary for the numerical computation due to the finite number of available energy 
levels. Thus N\ ong (E) is given by this additional term. For the exact formula we refer the 
reader to [22]. There are no proofs concerning the ergodic properties, but numerically the 
system behaves like a chaotic one. 

4. Isospectral billiard 

In 1966 Kac posed the famous question "Can one hear the shape of a drum?" [86]. Since 
Weyl's law determines the asymptotic behaviour of the spectrum in terms of geometric 
quantities like the area, circumference and properties of the corners of the membrane, 
the question arises whether there exist two billiards with identical spectra but different 
boundaries. Surprisingly there are indeed isospectral billiards as is shown in [87], where 
two planar isospectral billiards are constructed, each consisting of seven triangles with 
lengths (1,1, y/2). One of them is shown in figure 4. Thus Kac's question is answered in the 
negative. The first 25 frequencies are determined experimentally in [88]. Using a mode- 
matching method the first 598 quantal levels are computed in [89]. It is demonstrated that 



N el (E) = - 
N e2 (E) = - 
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the spectral statistics lies between the GOE and the Poisson behaviour with respect to the 
short-range correlations. However, when the quantal spectrum belonging to the triangle is 
subtracted, an agreement with GOE is observed. We do not pursue this approach because 
the spectrum of the triangle is not a subspectrum of the isospectral billiard since it cannot 
be obtained by desymmetrization. 

Using the boundary element method the quantal levels E n in the low-energy range E n £ 
[0,5000] and a high-energy range E n £ [37000,45000], have been computed. The first 
energy range comprises of 1335 quantal levels, the second one of 2210 levels. 

5. Triangular billiard 

Since the isospectral billiard is built up from triangles with lengths (1,1, a/2) , we want 
to study the eigenvalue statistics of one of these integrable triangles separately. The 
eigenfunctions are given by 

^m,n(x, y) = [sin(7rmx) sin(7rny) — s'm(iTnx) sin(7rmy)] 

where n,m£{l,2,...} with n < m. Thus the eigenvalues are 

E n ,m = 7T 2 (n 2 + ITT- 2 ) n, m £ {1, 2, . . .}, with n < m . 

The corresponding constants of the generalized Weyl's law in eq. (4) are given by A = |, 
C = 2 + y/2 and C = |. Altogether the first 10 6 eigenvalues are used in the following 
calculations. 

Note that for this system a non-Gaussian behaviour for the mode-fluctuation distribution 
P(W) holds, see section VI. 

6. Circular billiard (even symmetry) 

In case of the integrable unit circular billiard the (non-normalized) eigenfunctions are 
given in polar coordinates by 

± n (r, ip) = J m (k m , n r) e ±mv with m = 0, 1, 2, . . . , m = 1, 2, 3, . . . 

The boundary condition n (l,<p) = leads to the eigenvalues E m{a = k 2 m n as the squares 
of the positive zeros k m , a of the Bessel functions, i.e., J m (k m;a ) = 0. For m ^ the 
solutions are twofold degenerate; we will consider the case of the desymmetrized circular 
billiard with even symmetry with respect to the reflection at the diameter. The set of 
eigenvalues is given by {E mjfl } } and each eigenvalue occurs with multiplicity one. The 
corresponding constants in eq. (4) are given by A = ^, C = it — 2 and C = — Altogether 
the first 10 5 eigenvalues for even symmetry are used in the following calculations. 

7. Torus 

As a further integrable system we study a system which is defined on the two-dimensional 
torus (see [27, 50]). The eigenvalues of the "displaced" Laplace operator 

H = (iV + a) 2 = -(d 2 x + d 2 ) + a\ + o\ + %a x d x + 2ia 2 d y , 

with (x,y) £ [0,27r] 2 together with periodic boundary conditions read 

E ni ,n 2 = {ni - «i) 2 + (n 2 - a 2 ) 2 ; n 1 ,n 2 eZ , «i,a 2 £[0,l[ . 



25 



The eigenfunctions are given by 

iTr/ \ in\x inoy 

W(x,yj = e e y . 

In the following analysis we use the first 10 6 eigenvalues for the case a = (0.3437, 0.4304), 
which produces a bimodal distribution of PiW) [27, 50, 62]. 

For this class of systems it is proven in [50] that the corresponding limit distribution is 
non-Gaussian. The case with a\ = a 2 = was treated before in 



8. Rectangular billiard 

The last integrable system used in our comparison between chaotic and non-chaotic sys- 
tems is a rectangular billiard with ratio a/b = 5/7, and with Dirichlet boundary conditions. 
The (non-normalized) eigenfunctions are given by 

( Kmx\ ( irny\ 

^ m , n \ x i y) = sm sm ~~T~ ' m,n . (71) 



The corresponding eigenvalues are E m , n = it 2 ( + fr J • This system is chosen because 
of a recent experiment using such a microwave cavity [90]. In the experiment the first 657 
eigenvalues were measured and it is found that the mode-fluctuation distribution PiW) 
is very close to a Gaussian though the system is non-chaotic. Here we base our analysis 
on the first 10 6 eigenvalues, demonstrating clear deviations from the Gaussian, see below. 
Notice, that also for this class of systems a non-Gaussian behaviour of the mode-fluctuation 
distribution PiW) is proven, see section VI. 



VII. 3 Numerical comparison of chaotic and non-chaotic systems 

Let us now turn to the numerical comparison of the chaotic and non-chaotic systems which we 
have introduced above. In the numerical computations we have calculated a superposition of 
the statistics for different values of T, such that all available eigenvalues are taken into account. 
At first figure 5 shows the E(k } L) functions in the case of the cardioid billiard using the first 
10 000 quantal levels of the odd symmetry class in comparison with the GOE behaviour. The 
range L = ... 10 is displayed in figure 5a, and L = 50 ... 60 in figure 5b. One observes a 
good agreement with the GOE behaviour up to L ~ 3. Thereafter the width of the E(k } L) 
functions saturates, whereas the width of the GOE curves increases logarithmically as discussed 
in section III. 2. The increasing width corresponds to a decreasing maximal height of the (nearly 
perfect) Gaussians which is exemplified in figure 5b where the GOE curves are lying below the 
curves belonging to the cardioid billiard. This means that the spectrum of the cardioid billiard 
is at large correlation lengths more rigid than a GOE spectrum. The saturation effect is more 
clearly visible in figure 6 where the range L = ... 60 is shown. The variations in the widths 
and in turn in the amplitudes are determined by the number variance E 2 (L), see eq.(34), which 
fluctuates around the saturation plateau. 

In order to check the agreement of the E(k } L) functions with a Gaussian, we map the 
E(k, L) functions onto each other. To that end we apply the substitution £ = where c and 
a denote the mean value and the standard deviation of the E(k } L) functions, i. e., E n (k } £)d£ = 
E(k } L)dL = E(k } £cr + c)od£. The transformed E(k } L) functions are displayed in figure 7 
for the cardioid billiard in a), the stadion billiard in b), the truncated hyperbola billiard in 
c), the hyperbolic octagon in d), the hyperbolic triangle A in e) and the arithmetic hyperbolic 
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Figure 5: The E(k } L) functions are shown for the cardioid billiard using the first 10 000 odd 
quantal levels (full curves) in comparison with the GOE behaviour (dotted curves). The two 
ranges L = ... 10 and L = 50 ... 60 are displayed. 
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Figure 6: The E(k } L) functions are shown for the cardioid billiard using the first 10 000 odd 
quantal levels in the range L = ... 60. 

triangle B in f). In these figures the E(k, L) functions from k = 75 up to k = 125 are mapped 
onto each other. In all cases a Gaussian with variance a 2 = 1 lies well in the cloud of the 
E(k, L) functions, indicating the good agreement with a Gaussian. Nevertheless, there are 
system dependent differences as a comparison between the cardioid billiard with the hyperbolic 
triangle B shows. 

In figure 8 the same mapping is applied to five non-chaotic systems again for E(k, L) from 
k = 75 up to k = 125 in the case of the isospectral billiard, and from k = 5000 up to k = 5050 
for the other systems. The high k values are necessary because the latter four systems reach 
their non-universality regime much later than k = 125. Being in the non-universality regime is, 
however, required for a comparison with the mode-fluctuation distribution PiW). The results 
for the isospectral billiard are displayed in figure 8a and 8b for the low and a high energy 
range, respectively, as discussed in the previous section. It is seen that the distributions are not 
stationary since they change significantly from the low-energy to the high-energy range, which 
is in contrast to the case of chaotic systems where the Gaussian is approached very fast. Figure 
8c displays the E(k, L) functions for the Euclidean triangle billiard, figure 8d for the circular 
billiard, figure 8e for the torus and figure 8f for the rectangular billiard. 

The E(k, L) functions for non-chaotic systems are found to be neither Gaussian nor Pois- 
sonian. This is important because if they were Poissonian they would also tend to a Gaussian 
in the high k limit, and thus lead to a Gaussian behaviour of W(x) as we have seen in sec- 
tion IV. Only E(k, L) functions showing a non-Gaussian structure for L — > oo can generate a 
non-Gaussian PiW). 

A possible qualitative measure for the deviations from the normalized Gaussian behaviour 
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Figure 7: For six chaotic systems the E(k } L) functions are normalized as discussed in the text 
for k = 75 up to 125 such that they are mapped onto each other. The curves are in good 
agreement with the normalized Gaussian shown as white curve. 
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Figure 8: For five non-chaotic systems the E(k, L) functions are normalized as discussed in the 
text such that they are mapped onto each other. The A;-values are chosen out of the saturation 
range. Two distinct energy ranges are shown for the isospectral billiard with k = 75.. 125. For 
the other systems we have chosen k = 5000.. 5050. The curves show clear deviations from the 
normalized Gaussian shown as white curve. 



30 




is provided by the following L 9 -norm 

noT q (k) := ( r d£ 

We have tested different L 9 -norms and found that the norm for q = 1 shows the differences 
between chaotic and non-chaotic systems slightly clearer than higher L 9 -norms. In figure 9 the 
results for our systems are shown. To display all the nor 9 (A;)'s in a single figure, the A;- values 
have to be shifted according to k n = k — i.e., for the chaotic systems as well as for the 
isospectral billiard we subtract kj, = 75 and for the integrable systems kj, = 5000 from the k- 
values. Using the ordinate k n = k — the norms of the chaotic systems are plotted in figure 9 
as solid curves and of the non-chaotic systems as dashed curves. One observes a clear tendency 
towards larger values for the non-chaotic systems. However, there is the exception of the circular 
billiard whose norm is at some values of k n of the same order as for chaotic systems. On the 
other hand, at k n ~ 25, the L 1 -norm of the hyperbolic octagon and the triangle B display values 
ranging into the order of non-chaotic ones. But for the other non-chaotic systems the L 9 -norm 
yields a possible measure to discriminate between a Gaussian and a non-Gaussian behaviour. 

Another measure is provided by the moments and cumulants of r]L(x) which can be expressed 
in terms of the E(k } L) functions having the advantage that they catch different aspects of the 
distribution separately. To emphasize this point, figures 10 and 11 show the seven lowest 
cumulants of r]L(x) Ck(L) for the same systems as in figures 7 and 8. Here the cumulants have 
been smoothed in order to stress their mean behaviour. Because of the normalization, the second 
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Figure 10: The first seven smoothed cumulants Ck(L) are displayed for six chaotic billiards. 

cumulant C2 is identical to one. In the case of the chaotic systems all considered cumulants 
fluctuate around zero which is in contrast to the non-chaotic case where some cumulants Ck have 
large non- vanishing values revealing systematic deviations from the Gaussian behaviour. E.g., 
one observes for the triangular, the circular and the rectangular billiard large non-zero values 
for C4 and C&. In other cases more pronounced fluctuations occur than for chaotic systems. 
The cumulants of t]l(x) are computed from the moments which in turn are computed from 
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Figure 11: The first seven smoothed cumulants Ck(L) are displayed for five non-chaotic billiards. 
Two distinct energy ranges are shown for the isospectral billiard. 
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the E{k,L) functions using 

CO 

m 3 (L) = ¥ E{k,L) 

k=0 

The cumulants are then obtained by means of (48) 



(73) 



Ck{L) 



In M{n) 



fj,=0 



In J2 m 3{ L )-r 



(74) 



11=0 



where the series expansion of the logarithm yields the dependence of Ck(L) on rrij(L) with j < k 
[44]. From this one obtains (see eq.(47)) the cumulants of t]l(x), i.e., C\ = (ci — L)/y/c^ and 

Cfc = Ck/c^ 2 for A; > 1. 
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Figure 12: _D(x) is shown for the cardioid billiard using c = 1.1. The inset shows the same 
curve in comparison with that for the circular billiard where the increase is proportional to y/x. 

Before we finally turn to the PiW) distribution, we present in figure 12 D(x) for the cardioid 
billiard, using the first 10 000 levels and c = 1.1. In the inset additionally D(x) of the integrable 
circular billiard is shown. The first possesses a slowly logarithmically increasing D(x), whereas 
the latter displays an increase proportional to \fx. We have checked that the dependence on c 
is very weak. This provides us with a convenient method for computing the distribution PiW) 
in contrast to earlier methods. These earlier methods required, firstly, the determination of the 
spectral rigidity and its saturation value A O0 (E) in dependence on the energy E, and secondly, 
a fit to A O0 (E) in order to obtain its energy dependence [11, 15, 28]. This fit function yields 
an approximation to the variance D(x) of the mode fluctuations. Using the definition (52) for 
D(x), one obtains W(x) solely depending on the parameter c. From a numerical point of view, 
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skewness 


kurtosis 


cardioid billiard 


-0.018 


-0.14 


stadium billiard 


0.006 


-0.15 


truncated hyperbola billiard 


-0.036 


-0.12 


hyperbolic octagon 


0.015 


-0.18 


hyperbolic triangle A 


0.026 


-0.20 


hyperbolic triangle B 


0.010 


-0.06 


isospectral billiard (low energy) 


-0.172 


0.01 


isospectral billiard (high energy 


-0.366 


-0.03 


triangular billiard 


-0.318 


-0.04 


circular billiard 


0.233 


-0.12 


torus 


-0.062 


-0.79 


rectangular billiard 


-0.218 


-0.14 



Table 1: The skewness and the kurtosis of the distribution PiW) are listed. 

this parameter should be close to one since the larger the value of c the larger is the part of 
the spectrum being excluded from the analysis of PiW). However, a value too close to one 
implies an averaging over a statistically too small interval. The optimal value of c is thus a 
compromise depending on the number of available quantal levels. We have tested values in the 
range c = 1.1 ... 1.5 and observe only slight differences in the mode-fluctuation distribution 
P(W). In figures 13 and 14 the distribution PiW) is shown for the same systems as in figures 7 
and 8. Here c = 1.1 has been used, with the exception of the hyperbolic octagon where c = 1.5 
has been used. This is due to the observation that this system shows a c-dependence at c = 1.1, 
which can occur if not enough quantal levels are available. The difference in the mode-fluctuation 
distribution PiW) between the chaotic and non-chaotic systems is obvious. The chaotic systems 
display a very good agreement with the Gaussian behaviour, whereas significant deviations are 
observed in the case of the non-chaotic systems. 

The mean of the distribution PiW) is of the order 0.003 and the variances are accurate to 
1 ± 0.01 in most cases. In table 1 the skewness and the kurtosis are listed for the 12 mode- 
fluctuation distributions PiW). Relatively small values are observed for the chaotic systems. 
The modulus of skewness of the circular, the triangular and the rectangular billiard as well as 
of the isospectral billiard in the high energy range is larger than 0.2. The large value of the 
skewness in the case of the integrable circular billiard is also nicely reflected in figure 3 for 
^Vfluc(^), where the fluctuations are larger on the positive side than on the negative one. This 
leads to a suppression of the tail on the negative side in PiW) in comparison with the Gaussian. 
For the torus, where the skewness is of the order observed in chaotic systems, the modulus of 
the kurtosis is with 0.8 very large. Thus these two measures also display the deviations from a 
Gaussian distribution in the case of non-chaotic systems. 

A widely used test for the validity of a statistical hypothesis for discrete (cumulative) dis- 
tributions is the Kolmogorov-Smirnov test. Since the mode-fluctuation distribution PiW) is a 
continuous statistic, we apply this test to the distribution of the 8 n , eq.(6), analogously nor- 
malized as W(x) but using for the computation of D(x) only the values at the quantal levels. 
This can be viewed as a discrete version of the distribution PiW) since the fluctuations are 
only evaluated at energies corresponding to quantal levels. In table 2 the significance niveaus 
for a Gaussian distribution are shown where for the chaotic and for the isospectral systems all 
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Figure 13: The mode-fluctuation distribution PiW) is displayed for six chaotic billiards. 
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Figure 14: The mode-fluctuation distribution PiW) is displayed for five non-chaotic billiards. 
Two distinct energy ranges are shown for the isospectral billiard. 
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available quantal levels are used. For the remaining integrable systems only the first 10 000 
levels are used to employ for the comparison roughly the same order of quantal levels. It is 
clearly seen that the significance niveaus are much larger in the case of chaotic systems. Using 
for the integrable systems the first 10 6 levels leads to significance niveaus of 0%. However, using 
only the first 1000 levels yields significance niveaus larger than 10% for the circular and for the 
rectangular billiard. This shows that the deviations from the Gaussian are extremely significant 
for the integrable systems in figure 14 where 10 6 level have been used. This is in contrast to 
the case of chaotic systems where at most 10 4 levels have been used for the PiW) statistic in 
figure 13. 





significance niveau 


cardioid billiard 


48% 


stadium billiard 


77% 


truncated hyperbola billiard 


99% 


hyperbolic octagon 


91% 


hyperbolic triangle A 


63% 


hyperbolic triangle B 


77% 


isospectral billiard (low energy) 


12% 


isospectral billiard (high energy) 


10" 9 % 


triangular billiard 


10" 33 % 


circular billiard 


0.02% 


torus 


10" 14 % 


rectangular billiard 


10" 5 % 



Table 2: The significance niveaus of the distribution of the 8 n with respect to a Gaussian are 
given. 



VIII Summary 

In this paper we carry out a detailed comparison of chaotic and non-chaotic systems with 
respect to the statistical properties of their quantal level spectra. Emphasis is put on the 
E(k } L) functions and the mode-fluctuation distribution PiW). 

The RMT properties of the E(k, L) functions are discussed, as well as their behaviour for 
chaotic and non-chaotic systems. A crucial property is the Gaussian behaviour of the E(k, L) 
functions in the case of chaotic systems, which in turn leads to a Gaussian distribution of 
P(W). The L 1 -norm of the deviation of the E(k, L) functions from the Gaussian shows a clear 
tendency towards smaller values in the case of the six chaotic billiards, in contrast to the five 
non-chaotic billiards considered by us. In addition, the cumulants Ck(L) of t]l(x) reveal the 
same property since they fluctuate with small amplitudes around zero for the chaotic systems, 
whereas at least some cumulants display significant non-zero values or large fluctuations in the 
case of the non-chaotic systems. 

The distribution of t]l(x) is studied, which is known to be Gaussian in RMT and also for 
a Poissonian random process for L — > oo. Furthermore, a Gaussian behaviour of the E(k } L) 
functions implies a Gaussian behaviour of the distribution of t]l(x). 
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Finally the mode-fluctuation distribution PiW) is studied. There are rigorous proofs that 
P(W) is non-Gaussian for several integrable systems, leading to the suggestion that they provide 
a quantity for the distinction between chaotic and non-chaotic systems, because numerical com- 
putations show a Gaussian behaviour for chaotic systems. To confirm the conjecture that 
strongly chaotic systems with only unstable and isolated periodic orbits display a Gaussian 
P(W) in contrast to non-chaotic ones, we carry out a statistical analysis based on six chaotic 
and five non-chaotic billiard systems. The distribution of W(x) := Nf\ uc (x) / \J ' D(x) is numeri- 
cally shown to be in very good agreement with a Gaussian in the case of the six chaotic systems, 
whereas for the five non-chaotic systems clear deviations are found. 

A Gaussian limit distribution of PiW) for strongly chaotic systems is also satisfactory from 
the conceptual point of view. The often studied nearest-neighbour level-spacing distribution 
shows for integrable systems a Poissonian behaviour, in general, leading to the fact that in- 
tegrable or regular systems possess the most uncorrelated spectra. If one wants to have the 
degree of classical regularity correspond to the properties of a quantum mechanical statistic, a 
Poissonian behaviour for regular systems can be seen as surprising. In the case of the normal- 
ized mode-fluctuation distribution PiW), despite their rigid spectra, strongly chaotic systems 
display over long-range correlations the most uncorrelated spectra, which is expressed by the 
fact that a Gaussian distribution has the largest possible spectral entropy. 

Thus our results suggest that a clear fingerprint of quantum chaos is given by a Gaussian 
behaviour of the mode-fluctuation distribution PiW). 
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